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Qv , Abstract. There is a plethora of dark energy parameterizations that can fit current 

supernovae la data. However, this data is only sensitive to redshifts up to order 
one. In fact, many of these parameterizations break down at higher redshifts. In this 
paper we study the effect of dark energy models on the formation of dark halos. We 
select a couple of dark energy parameterizations which are sensible at high redshifts 
C^^ ■ and compute their effect on the evolution of density perturbations in the linear and 

^"T I non-linear regimes. Using the Press-Schechter formalism we show that they produce 

f~^ ' distinguishable signatures in the number counts of dark halos. Therefore, future 

VO i observations of galaxy clusters can provide complementary constraints on the behavior 

^^ ' of dark energy. 
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1. Introduction 



There are strong and converging evidences suggesting that the universe is nearly 
flat and that the dominant component of the energy density today has an unknown 
nature. Analysis of the luminosity distance of high redshift Type la supernovae has led 
to the conclusion that the expansion of the universe is accelerating, with indications 
that this acceleration is recent [1, 2, 3, 4, 5] and are corroborated by cosmic microwave 
background radiation [6, 7] and large scale structure observations [8, 9]. 

This suggests that the dominant contribution to the present-day energy budget can 
be described by a fluid with equation of state uj < —1/3, usually called "dark energy". 
We do not know for sure what this dark energy really is and in our ignorance many 
parameterizations for the time or redshift dependence of the equation of state uj{z) have 
been introduced in order to observationally constrain its evolution. 

Type la supernovae data are sensitive to uj{z) only for a small range of redshifts, 
typically up to 2; = 0{1). This is the reason for the large allowed regions of parameters 
related to the variations o{oo{z). In fact, the simplest possibility, a cosmological constant 
with oj = —1 still provides a good flt to the data. Therefore, it would be highly desirable 
to have another observable, sensitive to higher redshifts, that could break the degeneracy 
among the several different proposed parameterizations. 

Dark energy has a dramatic effect on the dynamics of the universe, altering the way 
in which cosmological structures grow [10]. This offers the possibility that observations 
of structure formation may provide a sensitive probe of dark energy properties that is 
complementary to information derived from supernovae data. 

The build-up of structure in universes with dark energy is subject to signiflcant 
ongoing study [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23], particularly through the 
spherical collapse model [24] in association with the Press-Schechter formalism for the 
mass function [25]. In this paper we use these techniques to investigate cosmologies in 
which the dark energy component remains homogeneous on the scales of the structures 
being modelled and their consequences to structure formation, in special to the cluster 
number counts. In particular, we show that cluster number counts can be sensitive to 
the different parameterizations of the dark energy equation of state. 

This paper is organized as follows. In section 2 we show the parameterizations we 
will be using to illustrate our results. In section 3 we study how the linear growth of 
perturbations is affected by these models of dark energy. Section 4 is devoted to the 
non-linear growth of perturbations in the spherical collapse model and the computation 
of the threshold density contrast for collapse in the different models studied. In section 
5 we present our results for the number counts of dark halos in different mass bins and 
also an integrated count in the different models studied. Finally, our conclusions are 
presented in section 6. 
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2. Parameterizations for the dark energy equation-of-state 



The expansion history of the universe is determined by the Hubble parameter, 
H{t) = a/ a, where a{t) is the scale factor (a(to) = aQ = 1 today). We assume that 
the dark energy has an equation of state relating its pressure pde and density poE at 
a particular instant determined by a{t) given by poE = ^{(i)Pde and that matter is 
pressureless, pm = 0. For general uj{a), the expansion rate of the universe is governed 
by the Friedman equation 

H^ia 






<) a-3 + 1]^ a-2 + fig;, e^(") , (1) 



where fi^, ^k and VtjJ^ are the current density parameters due to non-relativistic 
matter (baryonic and non-baryonic) , curvature and dark energy. In the following we 
assume a flat universe, VLk = (1 — ^m ~ ^de) = 0. Hq is the Hubble constant and the 
function /(a) is determined by the dark energy equation of state: 



and 



/(a) = 3 / [I + uo{a')]d\Yia' . (2) 

J a 

The matter density ^mI^) and dark energy density fiDB(o) are functions of a: 

fiM(a) = <)a-^-^, (3) 

^DE{a) = n%ef^^^^ (4) 



Throughout this paper we will adopt fiy = 0.25, fi]^;^ = 0.75, Hq = 72 km s~^ Mpc^"^ 
[7] and use either the scale parameter a or the corresponding redshift z = {1 — a) /a to 
describe the evolution of the different parameters. 

Recently discovered Type la supernovae (SNela) [4, 5] provide conclusive evidence of 
the decelerating universe in the past {z > 0.5) evolving into the present day accelerating 
universe. Thus the existence of dark energy, which accelerates the cosmic expansion, 
has been firmly established and the magnitude of its energy density today has been 
accurately measured. The goal is now to determine the behavior of the dark energy 
density and its equation of state at different cosmic epochs. The simplest and most 
natural candidate for dark energy is the cosmological constant A with a constant energy 
density ^a and a fixed equation of state parameter u = —1. It still provides a good fit 
to SNela data. 

The exact functional form of the equation of state uj{z) should follow from a 
fundamental theory for the dark energy. The best studied case is the so-called 
quintessence models, where dark energy results from a scalar field rolling down a 
potential. In this class of models, the equation of state is strictly in the region 
— 1 < a; < 1. However, more complicated models with multiple scalar fields or non- 
canonical kinetic energy can have uj < —1. 
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In the absence of a physically well-motivated fundamental theory for dark energy, 
it has become common practice to adopt parametric forms of uj{z) and to use SNela 
data to find the allowed regions in the parameter space for different parameterizations 
[15, 26, 27, 28, 29, 30, 31]. 

Many parameterizations were proposed to fit the observations and some of them 
are shown in Table 1. 

Table 1. Different parameterizations for dark energy (some of these models 
are from [31]) 

Model H^(z) or uj(z) Parameters 

la uj{z) = ujo + uJijY^^ = ujQ + uji{l - a)a c^o = -1-3 

c^i = 4 
lb uj{z) = ujQ + ujijj:^ = ujQ + uji{l - a)a c^o = -1-3 

ui = -2 

II u{z) =uJo + uJi (^^^)i.8 = cJo + cJi(l - a)a^-^ ujq = -IS 

uji = 4 

III H\z) = Hi [<)(1 + zf + n%\+ a, = 0.13 
+ai(l + 2;)^[cos(a2-2 + asTr) — cos(a37r)]] 02 = 6.83 

as = 4.57 

IV H^z) = Hl{n^^\l + zf + ai{l + z)+ ai = -4.16 
02(1 + ^)2 + (1^2)^-01 -as)} 02 = 1.67 

V H^z) = //o'[^M (1 + z?- Vai + 02(1 + ^)3 ai = 29.08 
(fig^ + Vai + «2)] 02 = -0.097 

VI uj{z) = ^^^^1^ ^0 = -1 

b = 0.25 

VII Uj{z) = UJq + Ui Z = Uq + Ui^—^ Uq = —1.4 

ui = 1.67 
VIII Lj{z) = uq + ^ij^ = ^0 + ^1(1 — cl) ^0 = —1.6 

ui = 3.3 
A u = -l or H\z) = Hl[n^^\l + zf + n^^'^E] - 

The parameterizations or models shown in Table 1 are obtained from several 
sources. Models la and lb are the same parameterization with the same central value of 
Wq but with values for wi at the extrema of the range allowed by SNLS and WMAP [32]. 
Model II is a slight modification of model I. Models III, IV and V are an oscillating, a 
quadratic polynomial [30] and a brane-motivated ansatze for H{z), respectively and we 
used the central values of the parameters [31]. Model VI was proposed by Wetterich 
[33] and its best fit parameters are taken from [18]. The parameterizations VII [26, 29] 
and VIII [15, 28] are first order Taylor expansions around z = and around a = 1, 
respectively. In figure 1 we show the behavior of the dark energy densities resulting from 
these parameterizations up to z = 5 and we see that most of them present problems 
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Figure 1. Evolution of the dark energy VLoe with redshift for parameterizations of 
Table 1. The lines are as follows: ACDM (solid), la (short dashed), lb (dotted), II (double 
dot dashed ), III (long dashed), IV (dot short dashed), V (dot long dashed), VI (double 
dot long dashed), VII (double dashed doted), VIM (triple dot dashed). 

i^DE(z) 




Figure 2. Evolution of the dark energy VLde with redshift for selected models of Table 1. 
The lines are as follows; ACDM (solid), la (short dashed), lb (dotted) and IV (dot short 
dashed). 



such as a too large dark energy component for redshifts higher than z = 1. 

For our study we select only a couple of models with decreasing contribution of 
^de{z) in the past. In order to contrast with the ACDM model, we choose models with 
energy densities both above and below the one given by a cosmological constant. These 
are models la, lb and IV, whose energy densities and equations of state are depicted in 
figure 2 and figure 3. 
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Figure 3. Equation of state evolution ll){z) with redshift for selected models of Table 1. 
Lines are the same as in figure 2. 



In the next sections we investigate the consequences of these different models to 
cosniological parameters that are relevant for large scale structure formation in the 
universe. 



3. Linear Perturbation Theory 



We assume that the dark energy component is smooth on scales smaller than the 
horizon [34]. In this case we only need to consider perturbations to the non-relativistic 
matter component. Dark energy only alters the background evolution and the equation 
of perturbed density contrast (for perturbations larger than the Jeans length) is [35] 



6 + 2H{t) 6 - -H{tfQM{t) 6 = Q, 



(5) 



where 5 is the fractional matter density perturbation. 

The growth function is defined as the ratio of the perturbation amplitude at some 
scale factor relative to some fixed scale factor, D = 6{a)/6{aQ), and its evolution equation 
can be written as 



D" + - 
2 



uj{a) 



X(a) 



D' 

a 



3 Xja) D_^ 
2 1 + X(a) a2 



(6) 



where the prime denotes derivative with respect to the scale factor a and X{a) is the 
ratio of the matter density to the dark energy density: 



X(a) 



n 



(0) 

M -3f^ dlna'uj{a') 



n 



(0)' 
M 



(7) 



For large X one recovers the matter dominated behavior D ^ a. The ratio of the 
matter density to the dark energy density X{z) can be seen in figure 4. Notice that for 
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all models except model IV the dark energy contribution decreases rapidly with redshift. 
In model IV the dark energy density decreases slowly but for redshifts larger than 0.5 
this parameterization does not accelerate the universe. 




12 3 4 5 

Figure 4. Behavior of the X with redshift for dark energy scenarios considered. Lines are 
the same as in figure 2. 

In figure 5 we show the evolution of the growth function D{z) for the selected 
dark energy models and for pure dark matter case. Notice that, as expected, larger 
perturbations in the past are needed to arise at the same amphtude today for models 
where dark energy is more important since negative pressure tends to inhibit the growth 
of perturbations. 

4. Non-Linear Evolution of Perturbations 



R 



(8) 



In order to describe the non-linear evolution of the density perturbations we adopt 
the spherical collapse model [24] where the radius R{t) of a spherical homogeneous 
overdensity region obeys the curvature-independent Raychaudhuri equation: 

-hl%{w{a) + l/3)e/('^)/? - hlf,{l + A,)-l, 

where time is measured in units of 1/i/o and Aj is the initial overdensity in the sphere. 
We numerically solved this equation for an initial time tj where a(tj) = 10~^, with 
initial conditions chosen so that the sphere is initially in the Hubble fiow, -R(tj) = a(tj), 
R(ti) = a(ti). We find the values of Aj such that the collapse occurs today. The linear 
evolution of 6i = Aj until today results in the critical linear density contrast parameter 
Sc which is important in the Press-Schechter [25] formalism discussed below. In figure 6 
we show an example in the ACDM model of the linear and non-linear evolution for an 
initial overdensity of Aj = 10"*^'^''^, chosen so that the collapse occurs today. 
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Figure 5. Evolution of the growth factor D{z) with redshift for selected dark energy 
models(lines are the same as in figure (2)) and for pure dark matter case (in which case 
D{z) = a{z) = yij) (long dashed line). 
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Figure 6. Evolution of the linear (dotted line) and non-linear (solid line) density contrasts 
in ACDM for an initial value chosen so that the collapse occurs today. 



In an Einstein-de Sitter universe, an exact value of 5c = 1.686 is obtained [24] and 
we verified that this value is fairly independent of the background cosmology [16]. 

We have also computed the values of Sc{z) for different collapse redshifts. Our 
results, shown in figure 7, find no large numerical differences among the models, in 
agreement with [10, 36]. We use these values of 6c{z) to compute the halo abundances 
in the next section. 
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Figure 7. Threshold overdensity for different redshifts of collapse for the models 
considered. Lines are the same as in figure 2. 



5. Mass Function and Cluster Number Counts 



The most reliable method to study the cluster abundance in the universe is through 
numerical simulations. However, there is an analytical approximation, the Press- 
Schechter formalism [25], that has been shown to fairly reproduce the results of N- 
body simulations [37]. There are more recent and better approximations with extra 
free parameters [38] but as an initial step we will use the well-known Press-Schechter 
approach. 

The basic premisses of the Press-Schechter formalism is to assume that the fraction 
of mass in the universe contained in gravitationally bound systems with masses greater 
than M is given by the fraction of space where the linearly evolved density contrast 
exceeds a threshold be-, defined in the previous section, and that the density contrast is 
normally distributed with zero mean and variance o"^(M), the root-mean-squared value 
of the density contrast b at scales containing a mass M. Therefore, it is assumed that for 
a massive sphere to undergo gravitational collapse at a redshift z its linear overdensity 
should exceed a threshold SJ^z\ Notice that only linear quantities are used in this 
formalism. 

These assumptions lead to the well-known analytical formula for the comoving 
number density of collapsed halos of mass in the range M and M -|- dM at a given 
redshift z: 



dn 



2 PmO ^c 



d\na{M,z) 



exp 



^l{z) 



2a^{M,z) 



(9) 



dM \ IT M a{M, z) dM 

where pmo is the present matter mean density of the universe and 6c{z) is the linearly 
extrapolated density threshold above which structures collapse, i.e., 5c{z) = 6{z = Zcoi). 
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The quantity 

(j{M,z) = D{z)aM (10) 

is the hnear theory rms density fluctuation in spheres of comoving radius R containing 
the mass M . The smoothing scale R is often specifled by the mass within the volume 
defined by the window function at the present time, see e.g. [39]. In our analysis we 
use the fit given by [12] 

where Mg = 6 x 1Q^'^Vl\J h^^ Mq is the mass inside a sphere of radius R^ = Sh^^Mpc, 
and (Tg is the variance of the over-density field smoothed on a scale of size R^. The 
index 7 is a function of the mass scale and the shape parameter, T = Vtjjh e^^^^^^/^M 
{Vtij = 0.05 is the baryonic density parameter), of the matter power spectrum [12] 



o"M = 0-8 



7(M) = 

Denoting 7(M) = 

7(M) = 

we can rewrite (9) as 
dn 



(0.3r + 0.2 

dlna(M,z) 



-n'-(^ 



dM 



(0.3r + 0.2 



(notice that it is z independent), 
2. f M 



2-^2+ 3 '°^M, 



2 PmO 5c{z) 

TV M a{M,z] 



■7(M) exp 



5ci4' 



2a{M,z) 



(12) 



(13) 



(14) 



dM 

In our study we use F = 0.144 [7]. For a fixed a^ (power spectrum normalization) the 
predicted number density of dark matter halos given by the above formula is uniquely 
affected by the dark energy models through the ratio 6c{z)/D{z). In order to compare 
the different models, we will normalize to mass function to the same value today, that 
is, we will require 

4,m(^ = 0) 



0"8,M 



Sc,a{z = 0) 



0"8,A 



(15) 
0.76 [7]. We show in figure 



where the label M indicates a given model and we use o"8,a 
8 the resulting mass functions for the different models. 

The effect of dark energy on the number of dark matter halos is studied by 
computing two quantities. The first is the all sky number of halos per unit of redshift, 
in a given mass bin 



dz 



dn 



4n 



^-p dn dV 
M. , dM dzdVt 



(16) 



where the comoving volume element is given by 

dV/dzdn = r\z)/H{z), (17) 

where r{z) = L H~^{x)dx is the comoving distance. The redshift evolution of the 
comoving volume element dV/dzdfi for different models of dark energy is shown in 
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(z) 

dM Mpc^ 




Figure 8. Press-Schechter mass functions for the different models with the cg 
normalization. Lines are the same as in figure 2. 

dv(z) Mpc^ 




Figure 9. Evolution of the comoving volume element with redshift for the four dark 
energy scenarios considered in this paper. Lines are the same as in figure 2. 



figure 9. Note that the comoving volume element does not depend on the growth factor 
of the perturbation D{z), but only on the cosmological background. The comoving 
volume element is larger for more negative equation-of-state, since this implies larger 
acceleration, see figure 3. 

The second quantity we compute is the all sky integrated number counts above a 
given mass threshold, Minf, and up to redshift z [20]: 

^'^- '■' > '"'"^ - L "^ Q, [m^ ''''' ■ (^«) 

Our knowledge of both these quantities for galaxy clusters will improve enormously 
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Figure 10. Evolution of number counts in mass bins with redshift for objects with mass 
within the range 10" < M/{h^ 
Lines are the same as in figure 2. 



within the range 10" < M/{h^'^MQ) < lO^"*. Notice the normalization factor of -^ 
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Figure 11. Evolution of number counts in mass bins with redshift for objects with mass 
within the range 10^* < M/(/i"1Mq) < 10". Notice the normalization factor of ^. 
Lines are the same as in figure 2. 



with upcoming cluster surveys operating at different wavebands [40] . 

The modifications caused by a dark energy component on the number of dark 
matter halos are tested and confronted with a cosmological constant ACDM model. We 
examine the effects of the different equations of state on the number of dark matter 
halos in mass bins [Minf,Msup] illustrating different classes of cosmological structures, 
namely [lO^^^ 10^^], [10l^ 10^^] and [10l^ 10^^] in units of H-^Mq. 

The number counts in mass bins, A/"*^™ = dN/dz, obtained from (16), are shown 
in figures 10, 11 and 12. In each of these figures we plot in the left panel the actual 
number counts and in the right panel we show the difference with the fiducial ACDM 
model. Notice that the more massive structures are less abundant and form at later 
times, as it should be in the hierarchical model of structure formation. Also there is a 
slight difference of the peak redshift for structure formation in the different dark energy 
models considered. 
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Figure 12. Same as figure 11 for objects with mass within the range 10^^ < 
Mjihr'^MQ) < 10^'^. Note that here TV and z range different from figure 11. Lines 
are the same as in figure 2. 
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Figure 13. Evolution of the integrated number counts up to redshift z for objects with 
mass above IO^'^H^^Mq. Lines are the same as in figure 2. 



The difference among tlie models results from a competition between the different 
volume elements and the different growth functions. At reshifts below one the comoving 
volume element has the most important role in the integral of Eq. (16). Above this 
redshift the comoving volume element does not vary significantly and the growth 
function becomes the dominant source for the number counts. 

An important observable quantity is the integrated number of collapsed structures 
above a given mass, equation (18). We present results for the integrated number 
counts of structures with masses above 10^^h~^MQ, IO^^H'^Mq, and 10^^h~^MQ. 
These are displayed in the figures 13, 14 and 15, together with the difference with 
respect to the fiducial ACDM model. We cut-off the integration in equation (18) at 
^sup = IO^^H^^Mq. Notice that the integrated number has a plateau that reflects the 
epoch of structure formation for a given mass. In other words, there is no formation 
of structures with mass above 10^^h~^MQ, 10^'^h~^MQ, and IO^^H'^Mq for redshifts 
roughly above z = 2, 1 and 0.6, respectively. 
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Figure 14. Evolution of the integrated number counts up to redshift z for objects with 
mass above 10^^/i~^Mq. Lines are the same as in figure 2. 
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Figure 15. Evolution of the integrated number counts up to redshift z for objects with 
mass above IQ^^h^^M^. Lines are the same as in figure 2. 



We notice that observations with accuracy of the order of 10%, either in binned 
or the integrated number counts, will able to distinguish among the different models, 
providing important information on the nature of dark energy. 

6. Conclusions 



We know that universe is currently accelerating and several different mechanisms 
have been proposed to explain the data. In the absence of a definitive model, an 
approach widely used in the literature is to assume a given parameterization of the dark 
energy equation of state. SNela data provides constraints in these parameterizations 
but are sensitive to redshifts up to 0{1). 

On the other hand, dark energy also influences the way in which the large scale 
cosmological structures form. In particular, number counts of collapsed structures is 
an important tool to probe dark energy models. Since structure formation occurs 
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at redshifts higher than those probed by SNela, number counts can provide useful 
complementary information in order to constraint the parameters of the different 
proposed parameterizations of the dark energy equation of state. 

In this paper we showed the impact of the different parameterizations and values for 
the parameters on several factors affecting large scale structure formation. We chose to 
exemplify our analysis by comparing a standard ACDM model with 2 parameterizations, 
one of them with two set of parameters currently allowed by data. 

We use the spherical collapse model in conjunction with the Press-Schechter mass 
function to investigate the effect of dark energy on the linear and non-linear growth of 
density perturbations and on the number counts of collapsed structures in different mass 
and redshift ranges. The dominant effect on number counts seem to arise from the mass 
function, which is more important than the comoving volume factor. The corrections 
arising from the merging of clusters were shown to be small [14] and hence were not 
considered in this work. 

Our results show that number counts can be useful in constraining dark energy 
models. Besides showing that the unintegrated and integrated number counts by 
themselves are powerful measurements, the redshift of maximum structure formation 
could also be used to differentiate the models. Observations with accuracy of the order 
of 10% can be used to distinguish among the different models, providing important 
information on the nature of dark energy. 

On the observational side, there are new experiments planned to start taking data 
on the abundance of clusters in the near future. In particular, the South Pole Telescope 
survey, based on the detection of the Sunyaev-Zeldovich effect arising from inverse 
Compton scattering of background photons off the hot intra-cluster gas, is expected 
to find a large number of clusters and will be able to determine the number counts to 
a high accuracy [40]. Hopefully these future observations will be able to discriminate 
among different parameterizations of dark energy currently proposed. 
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